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It is siiown tiiat a minimum realization of the dynamical mean-field theory (DMFT) can be 
achieved by mapping a correlated lattice model onto an impurity model in which the impurity is 
coupled to an uncorrelated bath that consists of a single site only. The two-site impurity model can 
be solved exactly. The mapping is approximate. The self-consistency conditions are constructed in 
a way that the resulting "two-site DMFT" reduces to the previously discussed linearized DMFT 
for the Mott transition. It is demonstrated that a reasonable description of the mean-field physics 
is possible with a minimum computational effort. This qualifies the simple two-site DMFT for a 
systematic study of more complex lattice models which cannot be treated by the full DMFT in a 
feasible way. To show the strengths and limitations of the new approach, the single-band Hubbard 
model is investigated in detail. The predictions of the two-site DMFT are compared with results of 
the full DMFT. Internal consistency checks are performed which concern the Luttinger sum rule, 
other Fermi-liquid relations and thermodynamic consistency. 






I 

o 
o 



> 

O 

in 

o 



a 



a 
o 
o 



X 



I. INTRODUCTION 



The dynamical mean-field theory (DMFT)la'Bffl has 
become a well-established and valuable method to inves- 
tigate the physics of strongly correlated electrons on a 
lattice. Similar as the Weiss mean-field theory for clas- 
sical models of localized spins, the DMFT is exact in 
the non-trivial limit of infinite spatial dimensions.™ For 
finite dimensions it provides a thermodynamically consis- 
tent and non-perturbative mean-field approach in which 
the spatial correlations are neglected but the temporal 
degrees of freedom are treated correctly. The apnlination 
of the DMFT to the single-band Hubbard modelQErEl has 
uncovered a complex phenomenology which may be char- 
acterized by strongly renormalized Fermi-liquid behav- 
ior competing with the Mott insulating state and with 
different kiaAa— nf .snontaneous order such as collective 



magnetism. 

The DMFT actually consists in a prescription that 
maps the original lattice model onto an effective impu- 
rity model which describes a single correlated impurity 
orbital embedded in an uncorrelated bath of conduction- 
band states. This mapping is a self-consistent one, 
namely the bath parameters depend on the on-site lat- 
tice Green function. The impurity model is the cru- 
cial point in the DMFT since it poses a highly non- 
trivial many-body problem that must be solved repeat- 
edly. The different methods employed for an essen- 
tially exact soluti|Qn-,of the impurity model, the, quan- 
tum Monte-CarlopE3 the exact diagonalizaticpifflij and 
the numerical renormalization-group method,ll30 work 
well for the single-band Hubbard model but are com- 
putationally expensive. In practice this severely limits 
the applicability of the DMFT, particularly when one is 
concerned with multi-band models. More complex lat- 
tice models including two, three or more degenerate or 
non-degenerate d-like bands possibly hybridized with un- 
correlated s-p-like bands are interesting for obvious rea- 
sons. For example, multi-band models are required for 
a minimum theory of the physics of strongly correlated 



electrons in the transition-metal oxides. Because of the 
large parameter space and the complexity of the effective 
impurity problem, a detailed and systematic calculation 
of the phase diagrams of multi-band models covering the 
entire parameter space is far beyond the ability of present 
implementations of the DMFT. 

If one is interested in a comprehensive mean-field de- 
scription of complex lattice models but wishes to keep 
the essence of the DMFT, some compromise is inevitable. 
In principle, there are only two possibilities conceivable: 
First, one may refrain from a numerically exact solu- 
tion of the impurity problem and employ approximate 
treatments instead. These may be based on different 
limits where a small parameter is available. This idea 
has been pursued with weak- and strong-coupling per- 
turbatiaaaLapproaches such as the iterative pertuthation 
theorytliiao or the non-crossing approximation.ESEElO In 
fact, multi-band Hubbard- type models can be treated in 
this way with an acceptable computing time (see Refs. 
^,^ for examples). This route, however, shall not be 
followed up here. 

The present paper takes into consideration the only al- 
ternative left, namely to solve the impurity model exactly 
but to reduce the number of bath degrees of freedom to 
keep the calculations manageable. There are no problems 
to treat e. g. the single-impurity Anderson model (SIAM) 
with a small number of sites rig numerically exact. On 
the other hand, for any finite rig < oo the self-consistent 
mapping of the Hubbard model onto the SIAM is ap- 
proximate. The exact solution of the effective impurity 
model is thus achieved at the expense of an approximate 
self-consistency. 

As a function of rig the Hilbert-space dimension D 
of the impurity model increases exponentially. It is 
given by D = 2"^^ where M is the number of (two- 
fold spin-degenerate) one-particle orbitals. The self- 
consistent mapping within the DMFT at least requires 
M = r -\- r(ns — 1) = rus orbitals for the case of one im- 
purity site, ns~ 1 bath sites, and r correlated bands (see 
appendix). Consider, for example, the d-band of a 3d 
transition metal with a two-fold degenerate Cg-derived 



band and a three-fold degenerate i2g-derived band. In 
this case M = ^Ug orbitals have to be considered. Ac- 
cepting M — 10 as a typical value for the maximum 
number of orbitals that can be taken into account for a 
repeated solution of the impurity model within a reason- 
able computing time, limits Us to the smallest number 
that is reasonable: ng = 2, i. e. an impurity model with 
one impurity site and one bath site. 

The purpose of the present paper is to test whether 
or not a two-site DMFT could be a meaningful approach 
and to show up its strengths and limitations. There are 
mainly two tasks to be performed: First, it is necessary 
to specify the approximate self-consistent mapping. This 
means to find a sensible prescription how to fix the bath 
parameters of the impurity model - guided by the origi- 
nal self-consistency condition of the DMFT. Second, the 
resulting two-site DMFT has to be tested against the 
full DMFT. This can be done best for the single-band 
Hubbard model (in infinite dimensions) where numerous 
essentially exact results are available. Thus, except for a 
short discussion of the extension of the theory to multi- 
band systems, the present paper is exclusively concerned 
with the single-band model. Since one cannot expect 
a strongly simplified two-site approach to reproduce the 
known results quantitatively exact, the present study at- 
taches importance to main trends, qualitative correct- 
ness and internal consistency. Recall that even the full 
DMFT cannot be expected to give more than a qualita- 
tive (mean-field) description of the physics of transition- 
metal oxides. 

The two-site SIAM has been considered-within the con- 
text of the DMFT beforehand. LangeEfl discussed the 
two-site SIAM at half-filling to investigate renormalized 
versus unrenormalized perturbational approaches to the 
DMFT of the Mott transition. However, a self-consistent 
mapping of the Hubbard model onto thcr-two-site SIAM 
was not considered. Bulla and Potthoflt3 developed a 
linearized DMFT of the Mott transition. The Hub- 
bard model at half-filling and at the critical interaction 
strength U = Uc was self-consistently mapped onto the 
two-site SIAM resulting in a linear algebraic mean-field 
equation. The solution of the two-site SIAM is exact, 
the mapping is approximate. It was found that, com- 
pared with the full DMFT, the linearized DMFT gives 
fairly good analytical estimates for Uc on different lat- 
tices. The theory, however, is restricted to the critical 
point. Further extensions of the linearized DMFT for 
the Mott transition have been developed for a periodic 
Anderson model by Held and liyJlaES and for a two-band 
Hubbard model by Ono et alE3 Extensions of the litu 
earized DMFT to the Hubbard model for thin filmsE2l 
and semi-infinite latticescj have convincingly shown that 
the main trends in the geometry dependence of Uc can 
be predicted safely. Here, it will be shown that a two- 
site DMFT can be constructed which is not bound to the 
(Mott) critical point but is able to access the entire pa- 
rameter space and that reduces to the linearized DMFT 
at half-filling and U = Uc- 



The paper is organized as follows: The next section II 
gives a brief review of the DMFT and introduces the self- 
consistent two-site approach. Section III presents a va- 



riety of results for the single-band Hubbard model start- 
ing with the Mott transition. Hereafter, the Fermi-liquid 
phase off half-filling is addressed, and the Luttinger sum 
rule, other Fermi-liquid relations and thcrmodynamical 
consistency are discussed. A discussion of the two-site 
DMFT in relation to other methods and the conclusions 



are given in section IV. The generalization of the theory 
to multi-band models is presented in an appendix. 



II. TWO-SITE DYNAMICAL MEAN-FIELD 
THEORY 

To be definite the single-band Hubbard model on the 
Bethe lattice with infinite connectivity g i-^ oo is consid- 
ered: 
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The hopping is assumed to be non-zero between nearest 
neighbors i and j only, t = —tij > is the nearest- 
neighbor hopping integral. The on-site hopping to = ta 
is set to to = to fix the energy zero. Furthermore, U is 
the on-site Coulomb interaction, cj^ creates an electron 
at the site i with spin a =1, |, and riicr = cl^Cia- With the 
usual scaling of the hopping, t = t* / ^/q and t* = const, 
the model is non-trivial in the limit g i— > oo.Ll Setting 
t* = 1 fixes the energy scale for the present paper. For 
a paramagnetic, spatially homogeneous phase the on-site 
Green function G{uj) = {{cia ', clcr)) i^ i^ given by: 
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where /i is the chemical potential, and S(w) is the self 
enftrgv which is local (k independent) in the limit q i— ; 
oo.qE3 poix) denotes the free density of states: 
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for \x\ < 2t* 
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The bandwidth is W 
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A. Dynamical mean-field theory 

The DMFT essentially rests on the observation that 
the local self-energy is given by a (skeleton-diagram) 
functional T, = S[G] of the on-site Green function that is 
universal for a large class of models. Consider, in partic- 
ular, the single- impurity Anderson model (SIAM): 



Himp = y^ eddjdq + Ud\did]di 



X] <^kal^aka+ ^ Vfe(4a/c(T +h.c.) , (4) 

a,k=2 CT,fc=2 



which describes an impurity orbital dj.|0) with one- 
particle energy e^ and on-site interaction U that is cou- 
pled via the hybridization Vk to a bath of rig — 1 non- 
interacting orbitals aJ-^lO) with energies e^. The impurity 



Green function G; 



Gimp(w) — 



"fco-l 
ap(w) 



{{da]dl))^ is given by 
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+ fi- Sd- A{uj) - Simp(t^) 



(5) 



where A(cj) — X]fe=2 ^fe^/(^ + A^~^fc) i^ t^*^ hybridization 
function, and I]iinp(<^) the impurity self-energy. As usual, 
Ed = ^0 = 0- The important point is that the functional 
S is the same as for the Hubbard model. Simp = '5[G'i,np], 
because the same type of skeleton diagrams occur in the 
expansion of Sjmp. Choosing the bath parameters e^ and 
Vk such that 



A{uj) =UJ + fi-€d- ^hnpiu) 



1 



GH' 



i. e. such that the DMFT self-consistency condition, 

Gimp(c^) = G{uj) , 
is fulfilled, then at once 

Simp(w) = S(a;) . 



(6) 



(7) 



(8) 



Therewith, the original lattice problem is mapped onto 
the SIAM and can be solved by the following iterative 
procedure: Starting with a guess for the local self-energy, 
the on-site lattice Green function is calculated from Eq. 
(|). Via Eq. (|) or Eq. (@) the Green function and the 
self-energy define the hybridization function A{uj) and 
thus the parameters of the effective SIAM. Finally, the 
impurity problem is solved to get a new estimate for the 
self-energy. The cycles have to be repeated until self- 
consistency is achieved. 



B. Two-site SIAM 

The self-consistency condition (M) can be fulfilled rig- 
orously only for rig i— > oo, i. e. for a bath with an infinite 
number of degrees of freedom. This leads to the usual 
SIAM which represents an involved many-body prob- 
lem. To simplify the problem and to construct a two-site 
DMFT, the case Ug = 2 is considered here, i. e. an ef- 
fective SIAM that consists of one impurity site and one 
bath site only. This represents the most simple bath con- 
ceivable. 

For Hg — 2 the site index is fixed to the value k = 2 in 
Eq. (0). There are only two independent bath parameters 
left, the one-particle energy of the bath site Ec = efc=2 and 
the hybridization strength V = Vk=2- The hybridization 
function is a one-pole function 



A(w) = V'^/{uj + 



M 



=) 



(9) 



and the free {U — 0) impurity Green function is a two- 
pole function: 



Ga(-) = ^. 



2r \UJ -\- ^ — € 



^J■ 



(10) 



with e = (ed — ec)/2 and r = Ve^ + V^. The inter- 
acting impurity Green function Ginip(w) has four poles 
and the self-energy Eimp(i^) two poles in general. Closed 
analytical expressions can be derived for the symmetric 
model at half-filling (see Ref. |2J, for example). For the 
non-symmetric case the model can be solved straightfor- 
wardly by numerical means without any problems. 

For clarity, the theory will be developed for the param- 
agnetic phase of single-band Hubbard model. However, 
it is rather straightforward to consider in essentially the 
same way also different magnetic phases and/or more 
complicated models such as multi-band Hubbard-type 
models, for example. 



C. Self-consistency 

For the two-site DMFT the Eqs. (|) - @ and Eq. M) 
are retained. The original self-consistency condition ([^, 
however, must be reformulated. This means to find two 
physically motivated (self-consistency) conditions to fix 
the bath parameters Ec and V. 

Consider first the limit of high frequencies w i— > oo. 
The exact self-energy of the impurity problem (jj) can be 
expanded in powers of 1/lu: 



S(w) = Uud + 



U^nd{l - rid) 



to 



+ 0{l/u;'), (11) 



where n^ = rido- is the average occupancy of the impurity 
orbital: 



1 r° 

(4d,)- / ImGn^p{uj + iQ+)dij. (12) 



rid 



Inserting the expansion InW) into (Ef) , one finds the follow- 
ing high-frequency expansion of the on-site lattice Green 
function: 



Gico) = - 



1 tp- i^L + Urid 

' 9 

UJ LO- 

M^°^ + (to - fi)^ + 2{to - fi)Und + W^rid 



0(l/c.4) 



W 



(13) 



Here M2 = J2j=ii^'ij — J dxx^po{x) is the variance of 
the non-interacting density of states (||). Eq. ( p^ ) can be 
compared with the exact high-frequency expansion which 
is available by calculating the first-^ion-trivial moments 
of the interacting density of states.^ One finds that Eq. 
(O) in fact represents the exact expansion provided that 



TT-imp = 2nd can be identified with the filhng n = {uii) + 
{fiii) of the lattice model. It is therefore required that 

n-mip = n , (14) 



where the band filling is calculated via 



2 r° 

— / IniG'(w + iO+)dw, 



(15) 



Eq. (Q) is the first self-consistency condition. The high- 
frequency behavior of G(w) is important for the occur- 
rence and for the correct weights and centers of gravity of 
the two high-frequency Hubbard excitations in the spec- 
trum. With Eqs. (|l2|), (|lj) and (|l^, an integral form 
of the original self-consistency equation (Q) is fulfilled: 
/_^ dw Im G = /_^ dio Im Gi,„p. 

Consider now the low- frequency limit w i-^ 0. The 
exact self-energy of the impurity problem (H) can be ex- 
panded in powers of w, 

E(CJ) =0 + 6^ + 0(^^2) ^ (;Lg) 

with constants a and h. The definition z = 1/(1 — b) will 
be convenient, i. e.: 

1 



^ 1 - dj:{o)/dLu ■ ^^"^^ 

For a metal z has the meaning of the quasi-particle 
weight. Neglecting terms of the order w^, and inserting 
into Eq. (|) yields G(w) = G^™*^) (u) for smaU uj where 
C(coh) j-^-j jg -(-j^g coherent part of the on-site Green func- 
tion defined as: 



G(™'')(c 



dx 



Po{x) 



dx 



uj + H — X — a — buj 
Po{x) 



LU — z(x — fi + a) 



(18) 



On the other hand, the coherent part of the impurity 
Green function is: 

1 






LU + fi - ed 



1 



buj 



^, ,, . (19) 
UJ — z(ed — fi + a^ Ay^)) 

Comparing the /li^/i-frequcncy expansions of the respec- 
tive coherent Green functions, 

z'^{to - 11 + a) 



G(=°h)(t^) = _ 
UJ 



z\Kq 



UJ'- 
'{to 



P 



and 



4;t^(^) 



+ Oil/iu*) 



LO LJ^ 

z'^V'^ + z^{ed - H + af 



(20) 



+ 0{1/lo^) , 



(21) 



leads to the second self-consistency condition: 



V' 



zM. 



(0) 



(22) 



Thereby, the original self-consistency equation (|^) is also 
fulfilled at low frequencies up to 0{uj) in an integral way, 
namely by referring to the weight, the center of gravity 
and the variance of the coherent quasi-particle peak. The 
equations (|l^ and ( p2[ ) reformulate the high- and the 
low-frequency range of the original self-consistency con- 
dition (^ in an integral, qualitative form and are thus 
well motivated. 



D. Calculations 

With the two conditions (n^ and (p2) the two bath pa- 
rameters are fixed and can be calculated self-consistently: 
Consider the model parameters fo, ti U, fj, and pa{x) to be 
given. Starting with a guess for Ec and V, the two-site 
impurity model (H) is well defined and can be solved to 
find the average occupancy of the impurity level nimp 
and the self-energy S((jj). The latter directly yields the 
quasi-particle weight z — {1 — d'E{0) / doj)~^ and from the 
condition (|22| ) a new value for the hybridization strength 
V. Inserting S(w) into Eq. (|^), gives the on-site lattice 
Green function G(w) and, via Eq. (|l5|), the filling n which 
has to be compared with ni,„p. Finally, a new value for Cc 
is chosen such that the difference n — riinip is reduced in 
the next cycle - according to condition (O). The cycles 
have to be iterated until ec and V are self-consistent (and 

In most cases, one finds a self-consistent set of bath 
parameters €c and V such that the ground state of the 
effective SIAM lies in the 6 dimensional subspace char- 
acterized by the total spin-dependent particle numbers 
N^ — Ni — I. To get the self-energy, it is convenient to 
calculate Giinp(w) from its Lehmann representation and 
to use the Dyson equation in reverse: 



E(c.) = cZi^y 



Giaip{uj) 



(23) 



where the free impurity Green function is taken from Eq. 
( [lO|) . The calculation of the eigenenergies and eigcnstates 
is straightforward, but even for rig = 2 this can only be 
done numerically in general. From the computational 
point of view, the most time-consuming step, however, 
is the calculation of the filling from Eqs. (g) and ([l^). 
Since the self-energy of the two-site impurity problem is 
a real two-pole function of the form 



S(aj) = ao + 



ai 



a2 



LU — UJl LU ~ UJ2 



(24) 



the filling can be calculated more directly by a single 
one-dimensional integration: 



n — 2 dLu po{uj + ft — T.{luJ) 



(25) 



p{u) = pq{u -\- ^ — S(w)) is the interacting density of 
states. Since even a repeated evaluation of this equation 
during the search for self-consistency is not very crucial, 
extremely fast numerical calculations can be performed. 

The two-site DMFT is obviously exact in the limits 
n = and n = 2. In the empty-band limit, for example, 
the impurity self-energy vanishes since n\ 



imp 



I ^ 0. 

Furthermore, it is exact in the band limit f/ = 0. In 
the atomic limit < = one finds ^ = (and arbitrary 
Ec) to be a self-consistent solution. Since n — riimp, the 
self-energy is given by 



Sh-iIcc-) = [/- + - 
2 uj 



U^{n/2){\~nl2) 
+ fi-to- [/(I - n/2) 



(26) 



which is the correct atomic-limit self-energy. Hence, the 
two-site DMFT is also exact for t = and, any filling. 
Since Eq. ( p6|) is the Hubbard-I self-energy,Q the theory 
reduces to the Hubbard-I approach whenever there is a 
self-consistent solution with V ~ 0. Actually, this is re- 
alized only at half-filling n = 1 and for sufhciently strong 
U (see below). 

Although the original lattice problem is mapped onto 
an effective impurity problem with a finite number of de- 
grees of freedom, the dependence of physical quantities, 
such as the filling n, on the model parameters will gener- 
ally be smooth. Consider, for example, the function n(fi) 
for a given U. The interacting impurity spectral function 
generally consists of 4 isolated 6 peaks of different weight. 
An infinitesimal change of fj, is unlikely to cause a finite 
jump of riimp since a change of nimp is in first place caused 
by a redistribution of spectral weight among the S peaks 
rather than by /x crossing a pole. Bearing in mind that 
the bath parameters itself depend on fi, the function n(p) 
can be continuous even in a large /i interval. In fact, it 
is found that the chemical potential never crosses a pole 
except for some extreme cases {n — 0, n ^ 2, V = 0). 



III. RESULTS 



The two-site DMFT provides a very simple, com- 
putationally fast and non-perturbative mean-field ap- 
proach to correlated lattice models. To show its 
strengths and limitations, numerical results will be pre- 
sented for the single-band Hubbard model which has 

pSMffliiiffiSfflQSESiSifil Thi^ win 

set the basis for a discussion of the approach and a com- 
parison with other methods in section IV. Multi-bands 



models are considered in the appendix. The calculations 
have been performed for the Bethe lattice with infinite 
connectivity q. The free density of states is given by Eq. 
(|). Its width is W = 4. 



A. Mott transition 

For the symmetric case of half-filling n — 1 and 
/i = U/2, the approach can be evaluated analytically. 
Particle-hole symmetry requires Ec = U/2 = p to en- 
sure the first self-consistency equation n = rii^p. Thus, 
only the hybridization strength V has to be calculated 
self-consistently. 

Let us first consider the critical point for the Mott 
transition, i. e. U t-^ Uc where the linearized DMFT is 
recovered (see Ref. |5|): Coming from the weak-coupling, 
metallic side, the critical interaction Uc is characterized 
by a vanishing quasi-particle weight z t—f which, ac- 
cording to Eq. (E3) implies V i-^ 0. In this limit two of 
the four poles of the impurity Green function are located 
near uj = zLU/2 while two poles lie close to w = 0. Re- 
ferring to the latter, the quasi-particle weight z can be 
calculated as the sum of their weights. A straightforward 
calculation to leading order in V yields: 



181^ 
C/2 



(27) 



On the other hand, another relation between z and V is 
given by the second self-consistency condition (^2|) . This 
implies 



V' 



M(0) ^ y2 



(28) 



namely a simple linear homogeneous mean-field equation 
as is characteristic for the linearized DMFT of Ref. |2^. 
A non-trivial solution of this equation is only possible for 
U = Uc where 



Uc = 6\ M, 



r(0) 



(29) 



The second moment of po{x) is easily calculated: Afj = 
^ if^. = qt'^ = r^ = 1, and therewith Uc = 6t* = 
1.5W. This result of the linearized DMFT is very close 
to the WWlt of the projective self-consistent method 
(PSCM)ya, Uc = 1A6W, and to the result of thejiumer- 
ical renormalization-group (NRG) calculationJlZl Uc = 
1A7W. The linearized DMFT is able to predict very 
reliably the main trends in the dependence of Uc on the 
lattice geometry rey| ci) | f ff r i systems with reduced transla- 
tional symmetryJHa'E3E§Ej 

The two-site DMFT is not restricted to the critical 
point U — Uc but more general. For p = U/2 and n = 1 
the self-eneijgy of the two-site SIAM can be calculated 
analytically:!^ 



n^) 



u 

'2 



f/2 



1 



1 



UJ-3V 



3V 



With Eq. (n^ this gives the quasi-particle weight: 

1 
^ ^ 1 + C/2/36t/2 ' 



(30) 



(31) 




FIG. 1: U dependence of the quasi-particle weight 

at half-filling as obtained wittiin the two-site DMFT, 
the Brinkman-Rice a«flp«aclS (dashed), the exact- 
diagonalization methocLflEj (pglid), and the numerical 
renormalization-group approachtZI (dots). 



consistent with Eq. ( P7|) ior V i-^ 0. Together with the 
self-consistency condition (22) one arrives at: 



V' 



M, 



(0) 



1 + C/2/36F2 



(32) 



This algebraic but non-linear mean-field equation has the 
self-consistent solution 



V 



- . nAo) 



t/2 



A'^"' - — 
36 



(33) 



r(0) 



for U < 6\/M^ ' = [/c and y = else. This yields the 
U dependence of the quasi-particle weight at half-filling: 

f/2 



Z= 1 



C/2 



(34) 



The resuljL-is the same as in the Gutzwiller variational 
approach :E3 



2BR - i - 772— 



U. 



(35) 



c,BR 



However, the Brinkman-Rice critical interaction f/c.BR = 
— 16/_ dx X po{x) « 6.79t* > 6t* is considerably 
stronger. 

The result of the two-site DMFT is compared with 
the Brinkman-Rice solution and with the results of nu- 
merical solutions of the|i|]il.DMFT by using the exact- 
diagtmalization mcthodtJ'cJ (ris = 8 sites) and the 
NRGIIZI in Fig. |l} There is a good agreement between 
the results of the ED and the NRG calculations except 
for interactions close to Uc where the energy scale of 
the quasi-particle resonance cannot be resolved reliably 
within the ED method. The two-site DMFT overesti- 
mates the quasi-particle weight in the whole U range. 



However, the qualitative agreement with the full DMFT 
is better than could be expected for the rather sim- 
ple approach and clearly improves on the result of the 
Gutzwiller method. 

The two-site DMFT interpolates between the trivial 
U = limit and U = Uc- For U — Uc it is a, reasonable 
approximation to neglect (i) the influence of the Hubbard 
bands at high frequencies on the low-frequency (quasi- 
particle) peak and (ii) the internal structure of the quasi- 
particle peak (see the discussion in Ref. ^^. These are 
just the basic assumptions of the linearized DMFT which 
the two-site approach reduces to for U = Uc- For < 
U < Uc these assumptions are less justified. Yet, the 
quadratic behavior of z{U) for C/ i-^ as well as the 
eventually linear behavior for U ^-^ Uc is consistent with 
the findings of the ED and NRG calculations. 

For U > Uc the self-consistent solution is given by 
V — which implies that the two-site DMFT reduces to 
the Hubbard-I approach in this case (Eq. (pq)). This is a 
crude description of the Mott insulator, even if compared 
with the Hubbard alloy-analogy solution and the itera- 
tive perturbation theory. The main deficiency is that the 
widths of the Hubbard bands are largely underestimated 
(see Ref. pa, for example). 



B. Fillings n^l 



For the symmetric case n = 1 the first self-consistency 



equation n 



''imp 



is fulfilled trivially by particle-hole 



symmetry since Cc — U/2 — fi, and only the hybridization 
strength V has to be determined self-consistently. For a 
thorough test of the two-site DMFT it is thus necessary 
to consider fillings off half-filling, too. 

Fig. shows the self-consistent bath parameters ec and 
y as a function of the filling for U — 4- As C/ is smaller 
than Uc the hybridization strength is finite for n — 1- 
For decreasing filling the system becomes less and less 
correlated, and consequently V has to increase until V = 
1 = M2 ior n — 0. According to Eq. ( |2^ ) this implies 
the correct value z = 1 for the quasi-particle weight in the 
empty-band limit. At half-filling the one-particle energy 
of the bath site is given by ec = 2 = U/2- It decreases 
with decreasing filling n and diverges on approaching the 
empty-band linnt, ec i-^ —00, as it is necessary to ensure 
a vanishing occupancy of the impurity orbital nimp = 
n i^ for finite V- Note that both parameters, ec and 
V, are smooth functions of the filling. 

Off half-filling the interacting impurity Green function 
continues to have four poles for U > and two poles for 
[7 = 0. This implies that the self-energy is a two-pole 
function of the form (|2J) not only for half-filling but also 
for n y^ 1- The poles wi and lu2 of the self-energy and 
their respective weights ai and 02 are shown in the lower 
panels of Fig. g. Again, it is noteworthy that these are 
smooth functions connecting the symmetric point n — 1 
where uji + uj2 = and ai = a2 with the empty-band 
limit where the poles become irrelevant since ai, a2 1-^ 0. 
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FIG. 2: Filling dependence of the parameters of the effective 
two-site impurity model ec and V and of the poles uji , uj? and 
respective weights ai, Q2 of the self-energy (see Eq. (P4)). 
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FIG. 3: Density of states p{uj) for U 
fillings. 



2W 



and different 



It is well knownQ that the density of states in the 
paramagnetic phase of the infinite-dimensional Hubbard 
model essentially consists of three peaks, the lower and 
the upper Hubbard band and a quasi-particle resonance 
in the vicinity of the Fermi energy. An attractive fea- 
ture of the two-site DMFT is that this general form of 
the density of states can be reproduced qualitatively. It 
is obvious that the two-pole structure of the self-energy 
results in a three-peak structure of the density of states 
(DOS) Pico) = poitu + fi - ^{uj)). 

Because the two-pole self-energy is real, the Hubbard 
bands and the quasi-particle resonance will be perfectly 
separated from each other on the real frequency axis. 
Clearly, this is only a sketch of the true density of states - 
the full DMFT generally predicts. the_i:espnance to merge 
with one of the Hubbard peaksD'EJ^OEj The symmetric 
case at half-filling is an exception. For U < Uc but close 
to C/c, a more or less clear separation of energy scales is 
found to be realized in fact.oO 

Results for the DOS are shown in Fig. g for U — 2W = 
8 at different fillings. Since U > Uc in this case, the Mott 
insulator is approached for n i~* 1. As soon as there is 
a finite hole concentration 1 — n, a quasi-particle reso- 
nance appears the width of which becomes broader with 
decreasing filling. The resonance is pinned to the Fermi 



energy. For the symmetric case n — \ and U < Uc it has 
a maximum at w = 0. For decreasing n < \ the maxi- 
mum shifts to a frequency wq > 0, and the asymmetry 
with respect to cj = increases. Concurrently, the up- 
per Hubbard band shifts to higher frequencies while its 
weight decreases. All this is qualit|atively correct when 
comparing with the full OMFTJaHEl However, the two- 
site DMFT largely underestimates the widths of the Hub- 
bard bands and does not predict the quasi-particle reso- 
nance to merge with the lower Hubbard band not even 
for smaller fillings. This is an obvious artifact. 

The width of the resonance is determined by z. The 
filling-dependence of the quasi-particle weight is shown in 
Fig. H for different U . For n i-^ there is a linear trend 
of the quasi-particle weight z — 1 oc n. For n i— j 1 the 
weight z{n) behaves linearly when U > Uc and quadrat- 
ically when U < Uc- Generally, z is a monotonously 
decreasing function of the filling for n < 1 and arbi- 
trary U and a monotonously decreasing function of U 
for arbitrary n which saturates in the limit U i—> oo. 
These resulis-are very similar to those of the Gutzwiller 
approachErEZrEa pjd qualitatively reproduce the results of 
the full DMFTE3 while quantitatively there are devia- 
tions similar as for the case n = I which has been dis- 
cussed already (see Fig. p. 
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FIG. 4: Quasi-particle weight as a function of the filling for 
different interaction strengths. Dashed line: calculation for 
[7 = 4 with a slightly modified theory (see text). 



The same holds for the fiUing dependence of the in- 
ternal energy E, the double occupancy d = {n^n\) and 
the chemical potential /i which is shown in Fig. |5|. For 
U = W = A as well as for U = 2W — 8 the chemical 
potential monotonously increases with increasing filling 
as it is required by thermodynamic stability. From the 
equation of motion for the on-site lattice Green function 
one can derive simple expressions for the kinetic energy 
per site, 



^kin — 2 



duj (lu + fi — 'E{uj))p{uj) 



and for the potential energy per site: 



E, 



pot 



div Y.{uj)p{u;) . 



(36) 



(37) 



The internal energy per site is given by _E = _Ekin + ^-pot, 
and the average number of doubly occupied sites by d = 
(n-^n^) — Epat/U . The double occupancy vanishes in the 
strong-coupling limit U i-^ oo. For U ~ W < Uc it is 
a monotonously increasing function of n while d is small 
(but non-zero) for U = 2W > Uc at half-filling. 



C. Test of Fermi-liquid relations 

Opposed to the fuU DMFT, the two-site DMFT 
is not a conserving approach in the sense of Baym 
and Kadanoff.E3 Consequently, the two-site DMFT can- 
not be expected to respect certain exact Fermi-liquid 
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FIG. 5: Filling dependence of the total energy E, the double 
occupancy d = (n^ni), and the chemical potential fi. Results 
for U = W (solid lines) and U ^ 2W (dashed lines). 



relations.^ Drastic violations of the Luttinger theorem 
and Fermi-liquid relations for the charge and spin sus- 
ceptibilities are well knowro for simple non-conserving 
theories such as the Hubbaui-I approximationJ3 the 
Roth two-ppple approximation,c3 or the spectral-density 
approach,c2l for example. For the two-site DMFT it is 
shown here that there are violations of Fermi-liquid re- 
lations indeed. Surprisingly, however, these are fairly 
small. 

One important property of a Fermi-liquid is that the 
self-energy T,{uj) is real at w = 0, i. e. fml](i0+) = 0. 
This guarantees the existence of a Fermi surface which 
(for a Bravais lattice) is defined by the set of points in 
k space with = ^ — e(k) — 2(0) where eji) is the free 
band dispersion. The Luttinger sum rulee^i states that 
the volume enclosed by the Fermi surface, 

yFS = 2^e(/i-e(k)-S(0)), (38) 

k 

is equal to the average particle number: 

^FS = (N) . (39) 




FIG. 6: Normalized "Fermi-surface volume" Vfs/ L as a func- 
tion of the filling n iox U — W (solid lines) and U — 2W 
(dashed lines). Results as obtained within the two-site DMFT 
and the Hubbard-I approximation. 



Here, 6 is the usual step function, and (N) — 
{J2ia '"'»o') = Ln where L is the number of lattice sites. 
The Fermi-surface volume can be calculated in the fol- 
lowing way: 



Vfs/L 



— oo 




dx po{x) Q{fJ. — X — S(0)) 

dx poix + fj. - J:iO)) (40) 



to be compared with the filling as given by Eg. (|25|) . 
The equation Vps/L = n with Vps/L given by (^) is a 
reformulation of the Luttinger sum rule that is also valid 
for the Bethe lattice. 

In Fig. o the Fermi-surface volume Vps/ L is compared 
with the filling n for different U. For moderate in- 
teractions ([/ = W, solid line) the Luttinger theorem 
Vps/L — n is almost exactly fulfilled in the whole range 
of fillings while deviations of a few percent are found for 
strong Coulomb interaction {U = 2W, dashed line) near 
half-filling. The results of the Hubbard-I approach are 
shown for comparison. It is seen that in this case the 
Fermi-surface volume is strongly overestimated up to a 
factor 2 near half- filling and irrespective of U. 

An alternative formulation of the Luttinger theorem is 
given byt2l 



/I = /IQ + S(0) 



where ^o is the chemical potential for [/ = 0. Replacing 
the first self-consistency equation (|lj) by Eq. (pl|), de- 
fines a variant of the two-site DMFT where the Luttinger 
theorem is enforced. The resulting filling dependence of 
the quasi-particle weight is shown in Fig. for [/ = 4 by 
the dashed line. If the original two-site DMFT respected 
the Luttinger theorem, there would not be any difference 
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FIG. 7: Charge susceptibility (compressibility) k and spin 
susceptibility x ^s functions of the filling for U — W = 4. 
Results for the two-site DMFT and the Hubbard-I approxi- 
mation as indicated. Solid lines: k and x calculated from the 
definitions (kd) and (|44[). Dashed lines: k and y calculated 
from the Fermi- liquid representations (|43|) and Wm. 



compared with the result of the variant. As is seen in 
Fig. U, the difference is non-zero but rather small. j—. 

By means of perturbation theory to all orderald the 
compressibility or charge susceptibility 



dn 
9/i 



(42) 



can be shown to be related to the DOS and the self-energy 
at ijj = 0: 



K = 2p(0) 1 






(41) Similarly, the spin susceptibility 



dm 
'db 



b=0 



is given by the expression 



X = P(0) 2- 



db 



(43) 



(44) 



(45) 
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FIG. 8: Filling dependence of the double occupancy d 
{n^ny) close to half-filling for different U . 



These equations provide two more exact Fermi-liquid re- 
lations. 

To calculate the spin susceptibiUty according to Eq. 
(E4), the formalism has to be generalized for the spin- 
polarized case. Consider a constant external magnetic 
field 6 in z direction which couples to the total spin: 



H ^ Ho- 2bS^ ^ Ho- bY^ 



qaTii. 



(46) 



where Ho = H{b = 0) is the Hamiltonian (|lj), riia- = 
cl^Cia, and qi — 1, qi — —1 is a. sign factor. The field 
strength is given by b. To account for a finite b in the 
formalism it is sufficient to redefine the on-site hopping 
to '^ toa- = to — Qa-b and likewise ed '-^ Cdo- — ^d — lab- 
Furthermore, the bath parameters may become spin- 
dependent: ec I— > £ctT and V t—f V^ ui the impurity model 
(^. The self-consistency conditions ([l4| ) and (^2|) now 
read: 



rida 



{dU.) 



vi 



,M. 



(0) 



(47) 



The magnetization is given hy m — n^^ — n<^. 

Fig. 7| shows the fiUing dependence of k and x as cal- 
culated from Eqs. (42) and (Kih t o be c omp ared with k 
and X as calculated from Eqs. (|43|) and ([45|). For U = A 
there are hardly any differences. Small differences of the 
order of a few percent are found for U = 8 (not shown) . 
On the other hand, within the Hubbard-I approach the 
filling-dependence of k and x strongly depends on the 
way it is calculated. The respective results are shown in 
Fig. 1^ for comparison. 

Concluding, one can state that the Luttinger theorem 
as well as certain Fermi-liquid relations are well respected 
by the two-site DMFT. The same holds for the ques- 
tion of thermodynamic consistency: Although the two- 
site DMFT cannot be derived from an explicit thermo- 



dynamic potential, consistency relations such as 



dUdn 



dU 



HI 



dn 



(48) 



are found to hold to a comparatively good approxima- 
tion. An example is shown in Fig. |g. At half-filling the 
chemical potential is given by /i = U/2. Therefore, as a 
function of the filling the double occupancy (n|nx) should 
have the slope d{n^ni)/dn = 1/2 at n = 1 for any (fixed) 
U. To a very good approximation this is reproduced by 
the results in Fig. ||. 



IV. CONCLUSIONS 

The two-site DMFT can be characterized as an approx- 
imate DMFT scheme which refers to a single-impurity 
Anderson model consisting of two sites only, one impu- 
rity and one bath site. Obviously, with a two-site SIAM 
the original self-consistency condition of the full DMFT 
can no longer be fulfilled, and consequently a compara- 
tively crude approximation has to be tolerated. The idea, 
however, is to construct in this way the most simple ap- 
proach that keeps the essence of the DMFT, namely the 
mapping onto an effective impurity model the bath pa- 
rameters of which are determined self-consistently. In 
fact, for the single-band Hubbard model some of the cal- 
culations can be performed even analytically and there 
are no serious problems to be expected for a numerical 
treatment of multi-band models. 

Any realization of the DMFT requires a repeated solu- 
tion of the impurity model which itself poses a complex 
many-body problem. An exact and unproblematic solu- 
tion is possible for a SIAM with a finite, small number 
of degrees of freedom - at the expense of an approximate 
mapping. This idea of the two-site DMFT is similar in 
spirit to the that of the exact-diagonalization method 
which aims at a minimization of the errors due to the 
discretization of the hybridization function by including 
as many sites Ug as feasible numerically. Contrary, the 
two-site DMFT stays with the case rig = 2 and can thus 
be considered as a "two-site ED" method. The latter, 
however, is not unique as it depends on the fit proce- 
dure used fioEj-the numerical determination of the bath 
parametersJljIij The advantage of the two-site DMFT is 
that it is based on well motivated self-consistency condi- 
tions to fix the bath. 

In principle, an extension of the two-site DMFT to- 
wards an rig-site DMFT is conceivable for the single- 
band Hubbard model. This requires the consideration 
of higher-order terms in the high- and low-frequency ex- 
pansion of the original self-consistency equation to fix 
the additional 2ns — 4 bath parameters. The compu- 
tational effort, however, grows exponentially with in- 
creasing Tig. Furthermore, higher-order static correla- 
tion functions would appear in the l/u expansion of 
the on-site lattice Green function (p^). At the order 
l/w^ the correlated-hopping correlation function B^ ex 
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Si5^j*u(cLciT(2ni_cr - 1)) can be expressed rigorously 
in terms of G{uj) by a sum ruleEil similar to Eq. (|l5|). At 
still higher orders, however, additional approximations 
must be tolerated. 

The implementation of the two-site DMFT is straight- 
forward, and the numerics is stable for the entire pa- 
rameter space. It has to be stressed that the approach 
allows for computations that are faster by several orders 
of magnitude compared to numerically exact approaches 
such as the QMC or ED method. This advantage quali- 
fies the approach for comprehensive investigations as are 
necessary e. g. for the determination of phase diagrams. 
In this respect it is coniO|a*able with the iterative pertur- 
bation theory (IPTjI^-Bta for example. 

The single-band Hubbard model has been considered 
here to illustrate the construction of the theory and to 
show up its advantages and limitations. Obviously, for 
the mere single-band model, the two-site DMFT cannot 
really compete against other more suitable methods. The 
main field of application are multi-band Hubbard-type 
models which require the solution of more complicated ef- 
fective impurity problems, or lattice models with reduced 
translational symmetry where oneriias-|to solve different 
impimty problems simultaneously I'^nrl Opposed to the 
IPTJ£3 for example, the extension of the two-site DMFT 
to the multi-band case is straightforward. Still the com- 
putational effort will be extremely small as compared to 
the numerically exact approaches. 

Taking the mass-enhancement factor m* /vi = z^^ 
as a measure for the strength of correlation effects, the 
Mott-transition point n = 1 and U = Uc obviously has 
a distinguished position in the phase diagram. Here 
the two-site DMFT reduces to the previously devel- 
oped linearized DMFTl3 which rather accurately pre- 
dicts (the meaDpfield) Uc for the Hubbard model on_dit. 
fercnt latticcsE^HEj and also for multi-band models E3'EZl 
The two-site DMFT can be considered as an extension 
of the linearized DMFT to the entire parameter space. 
It interpolates between the Mott-transition point where 
z — Q and the uncorrelated limit U = Q or n — Q where 
z = 1. Similar as the linearized theory, however, it can- 
not be controlled by a "small parameter" but is based on 
a physically motivated approximation. 

The two-site DMFT is the most simple approach that 
describes the transition from the Mott insulating state to 
the Fermi liquid as a bifurcation scenarios At half-filling 
and for U > Uc the approach reduces to the Hubbard-I 
approximation and yields the Mott-insulating solution. 
At U = Uc a, metallic solution splits off from the insu- 
lating one, the former being stable for U < Uc- Below 
Uc the two-site DMFT predicts a Fermi-liquid state with 
z ~ 1 — U'^/U^. TMa is the same as the Brinkman- 
Rice-Gutzwiller resultEj (albeit with a different Uc). Al- 
though conceptually the two-alaJDMFT is quite different 
from the Gutzwiller method,uO the results for depen- 
dence oi z on U and n are very similar. In particular, 
both approaches yield z{S) — z{0) ex S'^ for U < Uc and 
z{S) — z{0) ex S for U > Uc in the limit 5 = 1 — n i-^ 



(see Fig. ^ and Ref. ||). 

Similar as the IPT but opposed to the Gutzwiller 
method, the two-site DMFT does not respect the Lut- 
tinger sum rule for the invariance of the Fermi-surface 
volume since it is not _a. conserving theory in the sense 
of Baym and KadanofiEj and cannot be derived from an 
explicit expression for a thermodynamic potential. On 
the other hand, it is remarkable that the deviations from 
the Luttinger sum rule and also the deviations from dif- 
ferent exact Fermi-liquid and thermodynamical consis- 
tency relations are rather small. The comparison with 
the Hubbard-I approximation (which drastically violates 
these relations) suggests that this is due to the appear- 
ance of the quasi-particle peak in the spectral function. 
The Hubbard-I approximation can be considered as a 
local impurity approximation: The exact self-energy of 
the single-site (atomic) model is taken as an approxima- 
tion for the self-energy of the lattice mpdei. This ap- 
proach and similar but improved theorie^3c3 yield high- 
frequency excitations such as the Hubbard bands but fail 
to reproduce Fermi-liquid properties at low frequencies. 
To get the (low-frequency) quasi-particle peak in addi- 
tion and to restore the Fermi-liquid physics qualitatively, 
it is sufficient to couple merely a single (bath) degree 
of freedom to the single-site model. This also holds for 
the multi-band case. A fully conserving and consistent 
theory, however, can only be obtained by coupling to an 
infinite number of bath sites as in the full DMFT. 

Concluding, the two-site DMFT is a simple but non- 
perturbative mean-field approach to correlated lattice 
models. Compared to the full DMFT and to exact Fermi- 
liquid and thermodynamic relations, it yields satisfac- 
tory results for the Mott transition and the Fermi-liquid 
phase in the single-band Hubbard model with a mini- 
mum computational effort. The quality of the method, 
when applied to different physical problems, it not clear 
a priori and has to be examined. In particular, a com- 
parison between the two-site and the full DMFT concern- 
ing magnetic order in the single-band model as well as 
the question of finite temperatures are intended for the 
future. Primarily, future applications shall address the 
manifestly complex phase diagrams for spin, charge and 
orbital order in multi-band Hubbard-type models. Here 
the two-site DMFT may serve to give a quick and compre- 
hensive though rough overview of the main physics which 
can complement more thorough but selective studies. 
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APPENDIX A: MULTI-BAND MODELS 

The most general multi-band model with on-site inter- 
action reads: 



H 



E 



tiii2aia2 '^iiaia^i2a2C 



(Al) 

In most cases it is sufficient to take into account inter- 
action parameters that are labeled by two indices, the 
direct interaction Uaa' = t/oo'oo' and the exchange in- 
teractions Jaa' — Uaa' a' a and J'^^, — Uaaa'a'- Here 

a = 1, ...,r is an orbital index, and r denotes the num- 
ber of orbitals. Fourier transformation of the hopping 
Ui'aa' to k space yields the free Hamilton matrix ^^^'(k) 
the eigenvalues of which represent the free band struc- 
ture: em(k) where m = l,...,r is the band index. The 
orbitally resolved free density of states is given by: 

PoA^) ^lY. \<l>cm{k)\H{uj - e„,(k)) , (A2) 

km 

where (j)am (k) is the a-th component of the TTi-th eigen- 
state at each k, namely X^oo' 4'am(M)*taa'(M)4>a'm'0^) = 

^m,\^)^mm' ■ 

Assuming the self-energy to be local, it is easy to see 
from the usual diagram expansion that it must be di- 
agonal with respect to the orbital index: Tiaa'ai^) = 
Saa'^aa-i^)- Considcr a cubic lattice and let a refer to 
orbitals with an angular dependence given by the cu- 
bic harmonics. The lattice symmetries then require the 
on-site (i = i') elements of the lattice Green function 

Gij'oov(w) = {{ciacr; 4, ^,^))uj to bc diagonal with re- 
spect to a: 



(-'iiaa'a{l^) — Ooo' trocr ('^j 



(A3) 



Using the Dyson equation, the on-site Green function can 
generally be written as 



G. 



,„{lj) = 7E 



UJ + ^-t{k)-'Ecr{w) 



(A4) 



Here "Sa-i^^) is the (k- independent) diagonal r x r self- 
energy matrix and t(k) the k-dependent (non-diagonal) 
free Hamilton matrix with the elements taa' (k) . 

Within the DMFT the model ( Al ) can be mapped onto 
the following impurity model: 



Hu 



^da ^ct(T^Ck.(7 



E 

aa' 01020304 
aa.k=2 



Uaia2a4,a3 "oio-"o2ct' '^"3'^' "q 



Y^kaa 



E ^^° 



xaO'kaa -|-h.C.) , 



(A5) 



■,fe=2 



with e^a 



Gii„p,o<T(w) = {{daa;ct^aa))uj IS givCU by 



The impurity Green function 



Gimp, OCT ('J-') — 



1 



LJ + // - edo - ^aaii^) - Einip,oCT(w) 



(A6) 

The hybridization function Aq,o-(w) = X)fc ^kaali'^ + /^ ~ 
EfcoCT) and the impurity self-energy are diagonal with re- 
spect to a. 

The lattice self-energy can be derived by func- 
tional differentiation from the Luttinger-ffiard func- 
tional, TTiiiia'aa{ito) = 5$/(5Gji'QQ'cr (iLj) .El For any 
finite-dimensional lattice the DMFT consists in the as- 
sumption that the self-energy be local. This implies that 
$ depends on the on-site propagator Gaai^) only. Hence, 
the functional $ and thus the functional S = 5^/5G are 
the same for both, the lattice model (S — S[G]) and the 
impurity model (Simp = 5 [Gimp])- One can proceed as 
for the single-band case: If the bath parameters tkaa and 
VfcoCT are chosen such that 



AQ^(a;) =uj + 11- tda- Eimp,oCT(w) 



1 



Gqct(w) 
i. e. such that the DMFT self-consistency condition 



, (A7) 



Gil 



,(w) = Gacr{^) 



(A8) 



is fulfilled, then at once Simp,QCT(^) — '^acr{'^), and the 
usual self-consistent procedure can be set up. 

If different orbitals are equivalent due to lattice sym- 
metries as, for example, the three t2g orbitals in a cu- 
bic d-band system, one can make use of some simplifi- 
cations: The self-energy Y^aaii-jj) = ^^(w) and the DOS 
Po,o(w) = po(u}) are independent of the orbital index, and 
the on-site Green function is simply given by: 



Gacrij^) = Gcr{L0) 



dx 



Poix) 



UJ + ^1 — X — Ect(^) 



(A9) 



The lattice model can be mapped onto a simpler impurity 
model with ekaa = ^ka and Vkaa- = Vkcj ■ In this case the 
hybridization function and the impurity self-energy are 
a- independent, consistent with Eq. ( [AT] ). 

The two-site DMFT is constructed straightforwardly. 
The number of (two-fold spin-degenerate) one-particle 
orbitals in the impurity model is M = r-f r(ns — 1) = rrig. 
For Us — 2 and for the general case of r non-equivalent 
orbitals there are (for each spin direction) 2r bath pa- 
rameters to be determined, the one-particle energies €aa 
and the hybridization strengths Vaa for a = 1, ..., r. The 
comparison of the high-frequency expansions of Sq,o-(w) 
and Simp, OCT (w) to lowest order leads to a first set of r 
self-consistency conditions: 



"-d,, 



{daadao 



(ct 



(AlO) 
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with Uaa = (— l/"") /_oo I™ Gacr{^ + «0+) diO . 

In the low-frequency Hmit the self-energy is expanded 
as T.aai'^) — Oqct + (1 — -Zai)"^ + ^{lo'^) where Zaa is 
the orbital-dependent quasi-particle weight. Inserting 
t he self-en ergy expansion up to the linear order into Eqs. 
(A4) and (A6), yields the respective coherent parts of the 
on-site lattice and the impurity Green function. Analo- 
gous to the single-band case one gets: 

^mp,a.M - — + -2 

,2 T/2 I ^3 



^aa^aa ' ^aay^Aa f^ i Claa ) 



ur^ 



+ 0{1/lj^) 



(All) 



The expansion of the coherent on-site lattice Green func- 
tion, however, is slightly different: 

^(coh)i' N _ ^gg- ZgcrV'Oa ~ M + 0-aa) 



Zacr{t0a ~ M + Oqo-) 






+ Oil/to^). 



(A12) 



From the comparison one obtains the second set of r self- 
consistency conditions: 



aa ~ / , ^a'ljT^ii'aa' 



(A13) 



This includes an additional coupling of the d iffere nt or - 
bitals. Both self-consistency conditions (AlO) and (A13) 
reduce to Eq. (p7|) for the single-band case. 
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